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ABSTRACT 

t 

A  theoretical  investigation  is  developed  herein  in  order  to 
estimate  the  permanent  transverse  deflections  of  rectangular  plates 
when  subjected  to  large  dynamic  loads.     The    influence  of  finite  de- 
flections or  geometry  changes  and  of  strain- rate  sensitivity  is  retained 
in  the  analysis  but  elastic  effects  are  disregarded.     The  particular 
case  of  a  fully  clamped  plate     acted  on  by  a  uniformly  distributed 
pressure  pulse  is  studied  in  detail.     It  is  observed  that  good  agreement 
between  the  theoretical  predictions  and  experimental  results  has  been 
obtained  for  rectangular  plates  made  from  strain- rate  sensitive 
material. 

Thesis  Supervisor 
Professor  Norman  Jones 
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NOTATION 


B  Semi-width  of  plate 

D  Material  constant  defined  by  equation  (  2   ) 

H  Thickness  of  plate 

I  poZ 

V  poZ  /(  p-Hpc)^2 

L  Semi- length  of  plate 

6"oH 
Mo  Static  limiting  plastic  moment  = 

M'o  Dynamic  limiting  plastic  moment,   defined  by 

equations  (4a- d) 
No  Static  limiting  plastic  axial  force  =    6"oH 

N'o  Dynamic  limiting  plastic  axial  force,   defined  by 

equations  (4a- d) 
T,T  Durations  of  parts  of  the  response 

Vo  Initial  velocity 

Wm  Permanent  transverse  deflection 

Wf  Permanent  transverse  deflection  if   ^  H 

p  Material  constant  defined  by    equation  (2) 

pc  Magnitude  of  static  collapse  pressure 

pi,  po        External  pressures  per  unit  surface  area  of  a  plate 
po  Magnitude  of  dynamic  pressure   pulse 

t,t2,to      Time 

w  Transverse  deflection  of  the  midplane  of  a  plate 

x,  x',y      Co-ordinates  defined  in  figure  2 
z  Co-ordinate  which  is  perpendicular  to  the  mid- plane  of 

an  initially  flat  plate 

V>  JB 

L 

c  ij  Direct  strain 


U 


Po 
Pc 
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0m  Relative  angular  rotation  rate  across  a  line  hinge 

Kij  Curvature 

Mass  per  unit  area  of  plate 
£o  tan4> 

Non-dimensional  times 
6"o  Yield  stress 

~C  Duration  of  a  rectangular  pressure  pulse 


angle  defined  in  figure  2 

at 


(•) 

So  Strain  at  middle  axis  of  plate's  section 


INTRODUCTION 

In  reference  [1]   Jones    investigated  the  influence  of  finite 
deflections  on  the  behavior  of  arbitrarily  shaped  plates  and  beams 
which  are  subjected  to  large  dynamic  loads.     In  his  analysis   Jones 
disregarded  any  elastic  effects,   strain- hardening  and  rate  sensitivity. 
It  has  been  shown  that  elastic  effects  can  be  neglected  in  the  computa- 
tion of  final  deformation  if  the  kinetic  energy  input  is  appreciably 
larger  than  the  maximum  elastic  strain  energy  the  plate  can  absorb. 
With  respect  to  material  properties,   experiments  have  indicated  that 
strain  rate  sensitivity  is  by  far  the  most  important  property  neglected 
in  rigid- perfectly  plastic  methods.     Thus,  while  experiments  with  rate 
insensitive  materials,    such  as  Aluminium    A16061T6,   show  remarkable 
agreement  with  theoretical   results,   in  the  case  of  mild  steel,    a 
notoriously  rate  sensitive  material,    significant  disagreement  between 
theory  and  experiment  takes  place,   the  net  result  being  an  underestima- 
tion on  the  ability  of  a  plate  to  sustain  a  given  dynamic  load.     Cowper 
and  Symonds  \-^>°)  and  Symonds  and  Jones  l^l   have  developed  an 
approximate  theoretical  method  to  examine  the  influence  of  rate  sensi- 
tivity in  beams.     Their  results  show  drastic  improvement,   as  compared 

\2    51 
with  experimental  evidence,   over  the  non-rate  sensitive  methods.      * 

A  general  approximate  theoretical  procedure,   which  retains  both 

finite   deflections  and  strain  rate   sensitivity,    is  presented  herein  for 

the   dynamic  behavior  of  arbitrarily   shaped,   rigid- perfectly  plastic 

plates.     This  procedure  is  essentially  an  extension  of  the  one  developed 

by  Jones   "I  in  order  to  include  strain  rate  effects.      The   particular 

case  of  a   fully  clamped  rectangular  plate  subjected  to  a  uniformly 

distributed  dynamic  pressure  pulse  is  studied  in  more  detail. 
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APPROXIMATE  ANALYSIS 

In  reference  [1]  Jones  shows  that  if  a  plate  (or  beam)  is 
divided  into  a  number  of  rigid  regions  separated  by  V     straight  line 
hinges,   each  of  length  lm,     then  an  energy  balance  yields  the  following 
relation:  f 

\(p3-uw)wdA     =    ^__     J  (Nw-  M)  Qmdlm  (1) 

A  J  m  =  l      lm 

where  w  is  the  transverse  deflection  at  the  hinge. 

The  stresses  at  the  hinges  due  to  the  axial  force  N  and  bending 
moment   M   act  on  a  plane  which  is  parallel  to  the  hinge  and  transverse 
to  the  mid  plane  of  the  plate.      0m   is  the  relative  angular  rotation  rate 
across  the  hinge  and  is  supposed  to  be  infinitesimal.     The  quantity 
(Nw  -  M)  Qm  ,   on  the  right  hand  side  of    equation  (1),   represents  the 
internal  energy  dissipation  per  unit  length  of  a  hinge.     This  dissipation 
function  will  depend  on  the  boundary  conditions  of  a  plate  and  on  the 
yield  condition  used. 

For  a  strain  rate  sensitive  material  Symonds  and  Jones  in 
reference  [3]  represent  the  dependence  of  the  dynamic  lower  yield 
stress  S(j>)  on  the  plastic  strain  rate   £,      by  the  following: 

where    6"0     is  the  static  yield  stress,   and  p  and  D  are  material 
constants.     This  relation  is  totally  empirical  and  approximate,   pro- 
viding good  representation  of  test  results  up  to  strain  rates  in  the 
neighborhood  of   1,000   sec"*.     If  formula  (2)  is  used  to  derive  dynamic 

plastic  bending  moment  and  axial  force  for  rate  sensitive  behavior, 

[41 
assuming  the  Navier  hypothesis  valid,   Symonds  and  Jones  show 

that: 

(see  next  page) 
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M  4  z2 


1  -—779 (4a) 


Mo1  H2 


N  z 

n?-    =     2  IT  (4b) 


where    Mo'    =    [    1  +^7"    ( ^Ijp)1^  ]    Mo 


No'     =   [     l+f"3^)  %    ]    No 


and         z    ^  — 


M     a      i        H 

Mo        2p+1    sz 
-fo£  =     1 


,Vp 


wAere    JHo  _     [  1  4  (jo)  /f=j  _M0 
K=     [  i  +  (Lf}  N0 

CLV\X     Z      -^   H 

^  2. 

K  =    curvature  rate 

z     =   distance  between  middle  and  neutral  axes  of  the  plate's 
cross-  section. 

The  moment  Mo'    and  axial  force   No1  are  to  be  interpreted  as 
the  dynamic  equivalents  of  the  static  values    Mo,   No. 

When  equations    (4a)    and  (4b)    are  combined  they  give  the 
yield  condition: 

\Ml  .     /N\2  YIELD 


■  i-f^2  *^°  <5) 


Mo'  "VNo'/  CONDITION 
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This  condition  is  the  maximum  normal   stress  yield  criterion,   para- 
bolic in  shape,   as  is  seen  in  figure  1. 


N/No 


-i.0 

figure  1 

It  should  be  pointed  out  that  equation   (1)  has  been  derived  for 
arbitrarily  shaped  plates  and  since  it  is  an  energy  balance  equation,    it 
can  be  used  with  any  kind  of  edge  conditions,    provided  that   plastic 
hinges  are   allowed  to  form  around  the  boundary  of  the  plate,    if 
necessary,   as  well  as  the  interior. 

If  a  kinematically  admissible  collapse  mechanism  with  straight 
line   hinges  is  postulated,    then  equation   (1),    combined  with  the 
appropriate  relations  for   (Nw  -  M)  Om,   can  be  solved  to  give  an 
estimate  of  the  influence  of  finite- deflections  and  strain  rate  on  the 
deflection- time  history  of  rigid,    perfectly  plastic  plates  (or  beams) 
loaded  dynamically. 


DYNAMIC  BEHAVIOR  OF  A  FULLY  CLAMPED  RECTANGULAR  PLATE 


From  infinitesimal  plasticity  analyses  it  is  evident  that  the 
shape  of  the  displacement  field  of  a  beam  or  plate  for  small  dynamic 
loads  is  the  same  as  the  static  collapse  velocity  profile  l1»^J.     In 
addition,   from  experimental  results  I-    *     J  it  is  seen  that  the  permanent 
deformed  profiles  of  rectangular  plates  loaded  dynamically  are  similar 
to  the  shape  of  the  velocity  field  used  by  Wood  i6J  in  static  analysis. 

Consider  a  rigid,    perfectly  plastic  rectangular  plate  of  length 

2L  and  width  2B,   fully  clamped  around  its  outer  boundary,   as  shown  in 

figure  2.    The  plate  is  subject  to  a  uniformly  distributed  dynamic  load 

with  the  pressure- time  history  shown  in  figure  3. 
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This  external  loading  may  be  expressed  in  the  form: 
P3     =    Po  for    O^t  (r 


P3     =    0 


for  Z   <C  t     ^     T  +  X 


(6a) 


(6b) 


where    T      and         (T  +  T    )  are  the  durations  of  the  pressure  pulse  and 
the  plate  response,   respectively. 


2&: 


J  /  /  /   /    /    /     ////////  S  ///   J  /  J  ]/  ,_ 


figure  2 


0.0 


I  =  RT 


0.0 


figure  3 


Assuming  Wood's   displacement  field,    as  discussed  above,  we 


have: 


w    = 


for  region  I 
and: 


w    = 


Wi(B  tan<£       -x') 
B  tan4> 


Wi    (B  -  y) 
B 


(7a) 


(7b) 


for  region  II.     Regions  I  and  II  are  indicated  in  figure  2   and  i  =  1 

refers  to  the  deflections  during  the  time  imterval    0  ■%  t  ^.  T   ,     while 

i  =  2  refers  to    Z  ^  t     ^    T  +  Z    . 

In  order  to  use  equation  (1),    its  right  hand  side  is  examined 

first,   Thus:  ,  x 

(see  next  page) 
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/ 


) 


B 
cos^3 


(Nw-  M)Cmdlm  =4 


m  =  l     lm 


J 


(Nw-  M)Gm 


dx' 


sindp 


+ 


«J 


L-B  tan4>  L 

(Nw-  M)6mdx  +  4    J        (Nw  -  M)  €>m  dx  +  4  J  M) 

o  o  o 


B 

jo 


(8) 


The  evaluation  of  the  right  hand  side  of  equation  (8)  requires  the  use  of 
equations  (7a,   and  7b)  together  with  the   yield  condition  and  the 

relations: 

2 


Mo    = 


€o  H' 


No    =     6"oH 


Values  for  the  curvature  rate  K  ,   rotation  rate   ©m    and  hinge  width  to 
be  used  in  each  integral  are  given  in  table  (1). 


1       Integral 

2        Integral 

rd 
3        Integral 

4       Integral 

Om 

Wl 

2W1 

Wl 

Wl 

B  sin 

B 

B 

B  tan 

K 

Wl 

• 

Wl 

Wl 

Wl 

2B1  sin 

Bl 

Bl 

Bl  tan 

Hinge 

Width 

21 

21 

1 

i-H 

Table  (1) 
It  should  be  pointed  out  that  in  the  third  and  fourth  integrals  w    =  0. 

If  the  above  information  together  with  equations  (6a),    (7a)  and 
(7b)  are  substituted  into  equation  (1)  then  it  is  straightforward  to  show 

that:  ,  v 

(see  next  page) 
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Wl    +    (hl  Wl2    +   V  ^//p   +  h3Wl2    =   dl 


(9) 


where:    0  ^    t^Z   ,  W^H 


5         ,     h 


h. 


2        ,     d 


b3-l 


H-B' 


1  12  Mo 


(2-  b  tanft    ) 
(1+  b  cot<*>    ) 


b  = 


B 


,  1  1    +    2     (1  -  b  tan4>    ) 

2        3H2        [  (1  +  b  cot4>    ) 


Po 


Pc 


=  % 


Pc 


12  Mo        (  1  +  b  cot<4>    ) 


B 


2  (  3  -  b  tan<*>    ) 


k  A1    b 

b4    ' 


(1  +    b  cot4>    ) 


b5    ~  J 


3H     (  1  +  b  cot4>    ) 


A 


C2 

L 

B 

3C2 

L 

2  2 

C     sin  4=     +  C      x  C      cos  4> 


2  sin<±>    cos4> 


_     3C     tan4>  +  3C   -C     sin24>     -C     .    cos4> 

& i ^ O 

2  sin  4*    cos4> 
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r        2p         /      H       x1/ 

2       2p  +  1     v2  BID  p 


1       2p  +  1     v  4  BID  sin^ 


c      =-2e_     H 1/ 

3       2p  +  1        2  BID  tan<£   ;    p 


■x/  2 

and  tanc|?    =        y  3  +  b      -  b        to  ensure  smallest  upper  bound  to  the 

[5] 
actual  static  collapse  pressure. 

Equation  (9)  obviously  is  very  difficult  to  solve  analytically. 
Given  the  fact  that  the  whole  procedure  herein  is  an  approximate  one, 
and  aiming  at  providing  the  designers'with  relatively  simple  results,    a 
drastic  simplification  will  be  made  by  setting  p  =  1,   D  =  500  sec      ,   albeit 
the  proper  values  of  p  and  D  for  mild  steel  are  5  and  40  respectively. 
Thus,   the  highly  non- linear  equation  (2)  is  linearized.     Even  so, 
however,    it  is  hoped  that  the  whole  procedure  will  give  sufficient  insight 
into  the  behavior  of  strain  rate  sensitive  plates. 

Now  equation  (9)  may  be  recast  as: 

W     +  (h    W   2  +  h   )  W  +  h    W   2     =    d  (9a) 

I  1  J.  Z  O  1  1 

A  series  solution  to  equation  (9a)  employing  the  method  of  undetermined 

coefficients,   and  which  satisfies  the  initial  conditions  Wi  =  W±  -    0    at 

t  =  0,  yields  the  results: 

W      =   W 

•  • 

and  W      =    W 

at  the  end  of  the  first  stage  of  motion    (t  =  t    ),  where: 
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-  at2  u  u22u33,44 

W1     =      dir  [  1  -      h2Z      +   h2  r       -  h2   ^       +  h2    "C 


1 


12  60  360 


dlh3r4     +    W!     -     Y^        +     h2h3  dlT  l 


60  420  2520  126 


u    a    2  5 

hldl     T 

-         ]  (9b) 


84 


i,2        2           r,3        3  u44 

1                                  h2x         h       T             h      X  h2     T 

W      =      dr    fl-—       +— ^ -  — +  — 

1                1  L     L           2                  6                         24  120 


h3dlT  +      h3dlT  -      h2      T 


+ 


20  120  720 


,      h2h3dl    X  5    -       h!dx2    X5      j  (9c) 

36  24 


A  study  of  the  second  stage  of  motion  (  X  ^    u  ^   -*  "**^- ) 
proceeds  in  a  manner  similar  to  that  outlined  above  for  the  first  stage, 
but  with  equations  (9b)  and  (9c)    as  initial  conditions,   and  with     t\=  0. 

Thus  the  differential  equation  to  be  solved  now  is: 

w2  +  (hlw22  +h2)w2  +h3w22    .    d2  (10) 


where       d       = •         (     )     = 

a2  bl  l    '  dt 
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A  series  solution  to  equation  (10)  employing  the  method  of  undetermined 
coefficients  and  satisfying  equations  (9b)  and  (9c)  gives  the  maximum 
value  of  the  transverse  deflection  of  the  permanently  deformed  plate  as: 

• 

Wm  -  W2(T)  =Zl  +^_lZ  +2ilV  +a3x3P3  +dj/ 

H  H  H  H  H  H  H 


where:      CL      =      d2         -  J_  (h     +  h    W   2  )  W     -  J_    hQ  W   2 
1  2  2  l  l        2        3      l 


cc3  =-  _1_  [  a2  (h2  +  hl  wx2 )  +  (h3  +  hx  wx  )  wx  wx] 


a4  -  -  [2a2    W1  .  (3h1W1  +  h3)    +  3CC3(h2  +hx  Wx2  )    + 

+  (h3  +hj  wx)  wx2] 

T 
and  J3  =  t=-  is  the  smallest  non- negative  root  of  the  poly- 

nomial equation: 

4tt     tV  +  3a3^2+2a2V^l     =  °  (12) 


Equation  (12)  is  obtained  from  the  requirement  that  Wi  =  0  at  "W=  T 

where  t     refers  to  second  stage  alone. 

Equation  (11)  remains  valid  provided    Wm    ^  1,  because  when 

Wm     \  H   • 

— tz~  fi  1    the  dissipation  expressions    (Nw  -  M)  0m     take  on  new 

values  (  since  then  N  =  No'  and  M  =  0)  on  those  portions  of  the  hinge 

lines  which  have  transverse  deflections  greater  than  the  plate  thickness. 

In  this  case,   a  time  dependent  rectangular  boundary  travels  outwards 

from  the  central  line  hinge  toward  the  edge  of  the  plate.     This  boundary 

always  has  a  deflection  W  =  H    and  divides  the  plate  into  two  regions: 
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an  inner  zone  with    W  ^>  H     and  an  outer  zone  with    W  <\  H.     The  dis- 
sipation relations  (Nw  -  M)  0m     must  be  evaluated  in  the  inner  region 
with  N  =  No',   M  =  0    while  in  the  outer  region  the  evaluation  which 
was  used  earlier  remains  unchanged.     Figure  4  indicates  the  two 
regions:  |- 2L       *j 


2B 


'{//////'///'<'/<' 


2L 


SS//SS/S/    // 


* 


W/tt  <i 

t 


W/H  >1 


6 7^ 


4= 
7-rr 


77  //  ?  //■///  /  /  /  7-77777? 


figure  4 
If  it  is  assumed  that    W/H  exceeds  one  when    t    \~C  ,   then  the  first 
stage  of  motion  given  by  equations  (9a,  b,   and  c)    remains  unchanged, 
while  the  second  stage  remains  valid  until  t  =  ^  *  when  W/H  =  1 
where  tQ    is  the  response  time  of  the  second  stage  alone.     Thus,   the 


non-dimensional  time 
W  =  H,    i.e.  : 


f: 


_  t2 


may  be  found  from  equation  (11)  with 


«4r4j>24  +«343  +  a2Ij'22tw^rwrH  =  0  (13) 


The  smallest   non-negative  root  of  ^13)  gives  j>„  ,    i.e.   the  response  time 

for  the  second  stage. 

Now  all  subsequent  behavior  of  the  plate  will  have    W  //  H    at  the 

centrally  located  plastic  hinge.     It  is  therefore  necessary  to  consider  a 

third  stage  of  motion  which  is  governed  by  the  dissipation  relations 

employing  N  =  No1,   M  =  0    where  necessary,   as  discussed  earlier.     It 

is  quite  straightforward  to  show  that  in  this  case  the  following  equation 

is  obtained:  .  > 

(see  next  page) 
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1/ 
W     +  Z    W     +  Z2      +  Z    W_W     Vp    +  z  .     W3       P        =  0      (14) 
J  W  add  4W 

3  3 

Again  the  simplification    p  =  1    is  employed  so  that  (14)  becomes: 

W3W3   +    (Z3W32  +  Z4)W3    +Z1W32+Z2     =0  (14a) 

where    t2(    t    <T+t2    ,  P3     =0 

Zj  =JLi_  ,     z2  =JLs_   ,    z3  =  v±_  ,    z4  -  ^5 

TTj  TT1  1T1  TTj 


TT  -      =         fc^  2  -    jo)  f 

1  Pc  (3  -Jo)  ;       *°    ~    btan* 


^2 


2  (3  -   3  lo  +  5o2) 
H  (3  -      o) 


TT      =        Hb      (cot4>    +  tan<£>  ) 
3  (1  +  b  cot4>  ) 


.-1, 


^         roof      ,   f  (  2  sin<#    /p   , 

H  (1  +  b  cot         ) 

C0bH      [cot4>(tan4>  )  "%       +  tan£  ] 
5  ^ 


(  1  +  b  cot4>    ) 


A  series  solution  to  equation  (14a)  with  the  initial  conditions  W     -  H, 

•  ■ 

W     =  W  at      t       =  0,  where  t      is  refering  to  third  stage  alone,   gives: 


W3-    H  +  W2t3+a;t32  +  a;t33  (15a) 


20- 


where: 


a2B-2H     [Z1H2  +  Z2+(Z4,Z3H2)W2] 


<=-^tZimV2+Z3HW22+a2(W2+Z4  +  Z3H2)^ 


and  W3  =    W2  +2«  t3   +   SOL^2  (15b) 


when    t3    =  Tj      ,      Wg =  0  ,      W3    =    Wf 


Substituting  for  T     in  equation  (15b)  and  retaining  its  smallest 
non- negative  root.T^  is  obtained,   i.e.   the  response  time  of  the  third 
stage  alone.     The  value  of  the  final  and  permanent  displacement  is 
given  then  by  equation  (15a)  with  to  =  T-*. 


H 


1  +   ^2  Tl         +    #2  Ti2       +        OC3  TjJ 


(16) 


H  H  H 

where  W2     is  obtained  from  equation  (13)  by  differentiation: 

w2  =  w1  *  2a2r  ?2  +  sa3r  \2  +  4a4r  3j>3         (16a) 

An  alternative  approximate  solution  to  equation  (14a),   for  the  third  stage 

of  motion,   which  is  valid  for  any  value  of  p,  but  for  only  a  short  range 

of  values  of  the  aspect  ratio    b,   can  be  obtained  as  follows:  *-  "o    "7 


Call 


W3 

w      - 
4             H 

Z'1=HZ1     , 

Z'3 

=  HZ3 

z.     .*» 

^  2           H      ' 

Z'4 

Z4 

H 
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z\ 

Z' 

o(     = 

1 

o<     = 

2 

1 

Z'        ' 

2 

Z', 

3 

4 

Then  equation  (14a)  becomes 


llw     z'  ,    .        '     1/ 


W4  +  Z'3W4(CX1  +  W4     p)  +  4         (<X2+W4i/p)  =  0    (17) 

W4 

By  inspection  the  constants  o<  and  °^0'  which  both  depend  on  the 
aspect  ratio  b,  are  approximately  equal  for  0  ^ b  ^.0.3,  their 
ratio  varying  between   0.  98  and  1.16. 

Thus,   by  the  substitution  0<  ^  =c^o  =  °^  ,   equation  (17)  is  further 
transformed  to: 

W  Z' 

4  +   Z'    W      + 4_    =  0  (17a) 

*      1/  W 

((X+W    7p)  4 


Employing  the  change  of  variables 

d  W 
W     =  W      4  and  integrating: 

d  W„ 
4 


•* 


W41_2/P  d^[(W41"3/p  dW 


•     1-1/         *  1 

W4        ^   dW4    "   i 

L 

1-5 


'p   dW  -   -T7 
4         F   uvv4         c< 

r  t  W      "% 

1      /    •_  1-4/         •  1     ,•     1-5/         •  1    J  _±__  dW. )]  + 


4 


7  • 

3  2 

+        n  W/    +    Z'      log       W     +  C  =  0                                          (17b) 

2  4                4        6  e        4 


For  mild  steel,  when  p  =  5,   equation  (17b)  yields: 

(see  next  page) 
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7  34  g  33 

32  31 

+  <-5-<*"2--5-*"F)w//5    +  <-5-<*  "T-50<"3)i>      + 

Z  /  4  b  4 


,-4  ~6  Z*3  2  4 

+  (0<         _o<     )w„       +         _  W/    +Z'     log       W     +C=0 

4  2  4  4         e 


The  constant   C   is  evaluated  with  the  conditions: 


at  t  =  t  W     =  Z3_   =  1,  W     =    ^2_ 

Z  4  H  4  H 


Then  the  conditions 


Wf 


at  t  =  T        W.   =  0,      W.    =— — -  give: 

1  4  4  H  s 


Z'  W„  W  7       W 


2  H  4        ee  H     '  2  '     H 

+(_5_  __5_   o<       5  _2_      /5  _5_  ^        5    _  _5_  ^  _2_       /5 

M        9  '        H    '  v  8  3  V    H     ' 

-2  -—  W         7,  -—  -3        W       6/r. 

+  (_5_o<        .J-*        5       )(_2_)7/5    +  (_5_<*    5  ^    3 2_      fc 

V2  7  VH  V  6  '  V     H 


4  R  W2  Z\ 

+  (c<-4  -  c*~b)   ( ±—  )  +   — i-  (18) 

H  2 


The  final  permanent  deflection  is  given  by  the  root  of  transcendental 
equation  (18)  and  can  be  expected  to  be  within  the  interval: 

(see  next  page) 
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w 


/   wf    / 


1  + 


-7    W 
1    2  v   H  .4 


34     W„     9/5 


5.  <X     5     _2 

9  H 


+(- 


-—  -1       W       8/^  o  -—     W       7/^ 

>Ao<       5        _5_  *         )e_I)/5+(A^2.lcX       5     _2       /5 
8  3  H    '  v2  7  M   H  ; 


31 


W       6/  W 

5CK"3)  ( —  )  /5     +  (CX  '4  .  cK-6)  ( — 2_)  i 

6  H  H        J 


-i    1 


+  (_Lo<     5 
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DISCUSSION 

The  approximate  method  presented  in  this  pe^er  has  been  used  to 
predict  the  permanent  transverse  deflections    (Wf/H)   presented  in 
Figures  (5-9)  for  fully  clamped,   rate- sensitive  rectangular  plates 
which  are  subjected  to  uniformly  distributed  impulsive  loading,   i.e. 
when  pT    =  jA-Vo. 

It  is  evident  from  figure  (5)  and  tables  (2-6,   9)  that  the  estimates 
which  were  obtained  using  equations  (11)  or  (16)  and  rectangular  plates 
with  b  =  0.  593  agree  reasonably  well  with  the  corresponding  experimen- 
tal values. 

It  is  apparent  from  figures  (6-    8)  that,   for  a  given  value  of 
impulse,   the  permanent  transverse  deflection  of  a  plate  is  essentially 
independent  of  the  magnitude  of  t\  when  T\_  is  larger  than      50.     It  is  also 
noticed  that  for  increasing  values  of  ti^  the  deflections  slightly  decrease, 
whereas  in  the  non- rate- sensitive  analysis  of  Ref.   [1]  they  are  seen  to 
slightly  increase.     The  probable  explanation  for  this  difference  in  be- 
havior could  be  that  in  the  rate- sensitive  case  with  increasing  l\  the 
rate- sensitive  effects  become  more  pronounced  the  net  result  being  the 
(favorable)  reduction  in  Wf/H.     However,    since  for  large  values  of  T^ 
the  deflections  are  essentially  independent  of  t\,   as  discussed  above, 
for  impulsive  loading,  when  \—*-  oo    ,   this  difference  in  behavior 
becomes  insignificant. 

From  figures  (9,    10)  and  tables  (9-11)  it  is  observed  that  the 
value  of  the  permanent  deflection  is  also  affected  by  the  value  of  the 
hinge  width  and  the  value  of  the  material  property  D.     It  should  be  pointed 
out  here  that  the  value  of  D  is  merely  a  value  used  in  equation  (2)  to 
describe  the   6"  —  £    curve  for  the  material.     The  proper  value  of  D  for 
mild  steel  is   40  sec       when  p  =  5.     The  mathematical  difficulties  that 
led  to  the  linearization  of  equation  (2)  by  setting  p  =  1  naturally  have  an 
effect  on  the  value  of  D  as  well,   thus  affecting  the  final  results.   However, 
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it  is  evident  from  figure  (  )  that  for  values  of  D  greater  than  300  sec 

the  variation  in  Wf/H   is  rather  small  so  that  if  D  =*  300  sec       is  used 
with  the  present  analysis  results  good  enough  for  engineering  predic- 
tions can  be  obtained.     Similarly,   the  optimum  value  of  the  hinge  width, 
which  is  unknown  at  the  outset,   may  be  seen  from  figure  (    )  to  be  1  =  5H. 
(These  observations  are  made  for  the  particular  case  of  mild  steel). 

All  the  theoretical  results  presented  in  this  paper  were  obtained 
by  assuming  that  angular  changes  across  the  plastic  hinges  were  suf- 
ficiently small  to  permit  tan<$>     to  be  replaced  by  0    rad.     This  simpli- 
fication provides  a  good  approximation  when  Wf/H  is  not  too  large. 
Moreover,   only  the  few  first  terms  were  retained  in  the  series  solutions 
of  the  differential  equations  encountered  in  the  analysis.     It  would  be 
necessary  to  consider  additional  terms  for  large  Wf/H  ratios  and  for 
cases  in  which  the  dynamic  to  static  pressure  ratio  (\)  is  small.     The 
yield  condition  derived  in  ref  [4]  and  presented  here  as  equation  (2)  is 
an  approximation  to  the  exact  one,   though  the  effect  of  this  approximation 
is  not  considered  to  be  important.     The  linearization  of  equation  (2)  by 
setting  p  =  1  is  considered  to  be  of  some  importance  since  then  a 
highly  non- linear  curve  is  replaced  by  a  straight  line.     However,   given 
the  fact  that  the  most  important  factor  in  dynamic  loading  of  plates  is 
the  influence  of  finite  deflections,  which  are  ignored  in  bending- only 
methods,   this  linearization  seems  to  be  acceptable, at  least  tentatively, 
since  the  results  are  reasonably  close  to  experimental  values  and   at 
the  same  time  a  good  insight  into  the  rate- sensitive  behavior  is  gained. 

The  approximate  procedure  presented  in  this  paper  could  be 
used  to  examine  the  influence  of  large  transverse  dynamic  loads  on  the 
behavior  of  beams  and  plates  which  have  any  shape  and  support  condi- 
tions,  provided  the  appropriate  plastic  hinge  pattern  is  postulated. 

Finally  it  should  be  remarked  that  the  predictions  of  this  analysis 
are  thought  to  be  valid  when  the  external  dynamic  energy  applied  to  a 
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beam  or  plate  is  appreciably  larger  than  the  amount  of  energy  which 
could  be  absorbed  in  an  elastic  manner. 

CONCLUSIONS    -     RECOMMENDATIONS 


An  approximate  theoretical  investigation  is  developed  herein  in 
order  to  estimate  the  permanent  deflections  of  rectangular  plates  and 
beams  subjected  to  large  dynamic  loads.     The  influence  of  geometry 
changes  and  of  strain- rate  sensitivity  is  retained  in  the  analysis  but 
no  elastic  effects  are  considered.     The  particular  case  of  a  fully 
clamped  rectangular  plate  acted  on  by  a  uniformly  distributed  pressure 
pulse  is  studied  in  detail.     Tables  and  graphs  are  provided   for  the  case 
of  fully  clamped  plates  and  beams  made  from  mild  steel  and  loaded 
impulsively.     It  is  observed  that  good  agreement  between  the  theoretical 
predictions  and  experimental  results  has  been  obtained  for  a  rectangular 
plate  of  aspect  ratio  b  =  0.593.      Improved  theoretical  predictions  are 
expected  to  be  obtained  if  equation  (11)  is  solved   for  the  correct  value 
of  p.     Though  equation  (18)  is  an  exact  solution  to  equation  (17),    it  is 
valid  only  for  a  short  range  of  aspect  ratios  (0  -   0.  3)  and  even  in  that 
range  it  contains  some  approximations.     Moreover,   equation  (17) 
corresponds  to  the  third  stage  of  motion  only  so  that  the  exact  solution 
(18)  should  not  be  used  unless  equation  (11)  has  been  solved  correctly 
as  well.     Exact  or  numerical  solutions  to  equations  (11)  and  (16)  would 
be  a  useful  complement  to  this  paper. 
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Equation 

Experimental 

(ll)or(16) 

Ref  [2] 

H-in 

Vo- ft/sec 

I' 

Wf/H 

Wf/H 

.0643 

134.08 

323.05 

3.806 

3.682 

3.542 

.0644 

152.60 

417. 16 

4.326 

4.253 

4.  120 

.0638 

165.06 

497.30 

4.723 

4.684 

4.650 

.0638 

180.06 

591.78 

5.  152 

5.  123 

5.  166 

.0647 

233.00 

963. 55 

6.574 

6.458 

6.420 

.0635 

234.08 

1009.604 

6.  730 

6.643 

6.  730 

.0998 

80.73 

51.64 

1.522 

.771 

1.046 

.0985 

118.95 

115.08 

2.272 

1.515 

1.890 

.0984 

161.70 

213. 10 

3.092 

2.690 

2.755 

.0983 

177.90 

258.46 

3.405 

3.053 

3.330 

.0982 

202.60 

335.90 

3.881 

3.564 

3.760 

.0984 

231.  13 

435.38 

4.420 

4.096 

4.300 

.1728 

69.69 

11.79 

.  72  7 

.339 

.310 

.  1728 

88.85 

19.  16 

.927 

.457 

.  515 

.  1728 

130.90 

41.  59 

1.366 

.678 

1.022 

.1728 

153.36 

57.09 

1.600 

.780 

1.257 

.  1728 

165.70 

66.64 

1.  730 

.833 

1.420 

.  1728 

166.26 

67.09 

1.735 

.835 

1.411 

.  1728 

171. 10 

71.06 

1.785 

.855 

1.580 

.1728 

178.02 

76.92 

1.857 

.883 

1.715 

Table  (2  ). 

b  -  0.593, 

1  =  5H,     D  = 

=  300  sec 

,     Tt=  200. 
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Equation 

Experimental, 

(ll)or(16)  Ref  [2]. 

H-in 

Vo- ft/sec 

Wf/H 

Wf/H 

.0643 

134.08 

232. 05 

3.562 

3.542 

.0644 

152.60 

417. 16 

4.  102 

4.  12 

.0638 

165.06 

497.30 

4.  505 

4.650 

.0638 

180. 06 

591. 78 

4.  915 

5.  166 

.0635 

234.08 

1009.60 

6.318 

6.730 

.0998 

80.  73 

51.64 

.763 

1.046 

.0985 

118. 95 

115. 08 

1.431 

1.890 

.0984 

124.20 

125.72 

1.62 

1.940 

.0982 

177.90 

258.98 

2.  935 

3.330 

.0982 

202.60 

335. 90 

3.412 

3.  760 

.0983 

216. 70 

383.49 

3.662 

4.  135 

.0984 

231.  13 

435.38 

3.908 

4.300 

.1728 

69.69 

11.  79 

.  336 

.  310 

.  1728 

88.85 

19.  16 

.452 

.515 

.  1725 

130. 90 

41.59 

.670 

1.  022 

.  1729 

165. 70 

66.57 

.816 

1.411 

.  1727 

178.02 

77.01 

.864 

1.  715 

Table  (3  ). 

Theoretical  vs 

>.    experimental  results,    b  =  0.  593,    I 

-1 


1  =  2.5H,  t^=  200. 
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Equation 

Experimental 

(ll)or(16) 

Ref  [2] 

H-in 

Vo-ft/sec 

Wf/H 

Wf/H 

.0643 

134.08 

323.05 

3.517 

3.542 

.0644 

152.60 

417. 16 

4.04 

4.  12 

.0638 

165.06 

497.30 

4.443 

4.650 

.0638 

180.06 

591. 78 

4.845 

5.  166 

.0647 

233.00 

963. 55 

6.043 

6.420 

.0635 

234.08 

1009.60 

6.220 

6.  730 

.0998 

80.73 

51.64 

.762 

1.046 

.0985 

118. 95 

115.08 

1.428 

1.890 

.0984 

124.20 

125.72 

1.586 

1.940 

.0982 

177. 90 

258.98 

2.932 

3.330 

.0982 

202.60 

335. 90 

3.408 

3.760 

.0983 

216. 70 

383.49 

3.658 

4.  135 

.0984 

231.  13 

435.38 

3.903 

4.300 

.  1728 

69.69 

11.79 

.332 

.  310 

.  1728 

88.85 

19.  16 

.443 

.515 

.1725 

130.90 

41.74 

.647 

1.022 

.  1729 

165. 70 

66.57 

.783 

1.411 

.1727 

178. 02 

77.  01 

.828 

1.715 

Table  (4). 

b  =  0.  593, 

D  =  460  sec" 

*,  li  =  200 
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H-in 

Vo-ft/sec 

.  1728 

69.69 

11.50 

.  1728 

88.85 

18.69 

.  1728 

153.36 

55.69 

.  1729 

165. 70 

64.93 

.1728 

166.26 

65.45 

.  1727 

178.02 

75.  12 

Table  (5). 

Theoretical 

results 

V  Wf/H 

.367  .0  747 

.486  .1077 

.808  .2145 

.872  .2327 

.876  .2338 

.938  .2508 


-1 


1^  =  200. 
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Equation 

Experimental 

(ll)or(16) 

Pef  [2] 

H-in 

Vo-ft/sec 

Wf/H 

'Vf/H 

.0998 

80.  73 

51.64 

.765 

1.046 

.0985 

118. 95 

115.08 

1.459 

1.890 

.0984 

124.20 

125. 72 

1.618 

1.  940 

.0982 

177. 90 

258. 98 

2.976 

3.330 

.0983 

216. 70 

383.49 

3.  718 

4.  135 

.0982 

202.60 

335.90 

3.462 

3.  760 

.0984 

231. 13 

435.38 

3.969 

4.300 

.  1728 

69.69 

11.  79 

.333 

.310 

.  1728 

88.85 

19.  16 

.445 

.  515 

.  1725 

130.90 

41.74 

.653 

1.022 

.  1729 

165. 70 

66.57 

.  791 

1.411 

.  1727 

178.02 

77.01 

.837 

1.  715 

Table  (6  ). 

b  =  0.593, 

D  =  500  sec" 

1    <^=  ZOO 
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H-in 

Vo-ft/sec 

r 

Wf/H 

.0998 

80.73 

51.64 

1.769 

.  9520 

.0985 

118.95 

115.08 

2.642 

2.461 

.0984 

124.20 

125.72 

2.761 

2.607 

.0982 

117.90 

258.98 

3.963 

4.  138 

.0982 

202.60 

335.89 

4.513 

4.  758 

.0984 

231. 13 

435.38 

5.  139 

5.424 

.0983 

216.70 

383.50 

4.823 

5.091 

Table  (7). 

b  =  0.75,   D  = 

-  500  sec      , 

1  =  2.  5H,  i 

i  =  200. 

H-in 

Vo-ft/sec 

r 

.0998 

80.73 

51.64 

2.074 

.0985 

118.95 

115.08 

3.097 

.0984 

124.20 

125.72 

3.237 

.0982 

117. 90 

258.98 

4.646 

.0982 

202.60 

335.89 

5.290 

.0984 

231. 13 

435.38 

6.023 

Table  (8  ). 

b  =  1.0,     D  = 

:  500  sec       , 

1  =  2. 

Wf/H 

1.931 
3.605 
3.786 
5.  773 
6.617 
7.540 
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Equation 

Experimental 

1          H-in 

(ll)or(16) 

Ref  [2] 

D-sec 

Vo- ft/sec 

Wf/H 

Wf/H 

500 

0.643 

134.08 

3.698 

3.542 

450 

3.504 

400 

3.437 

350 

3.356 

300 

3.255 

250 

3.  128 

200 

2.96 

150 

2.728 

100 

2.376 

75 

2.  113 

60 

1.896 

40 

1.393 

Table  (9). 

Influer 

ice 

of  D  on  predict 

ions  based  on 

equations  (11) 

(16).     1  =  2.5H,  *i=  200. 
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1/H  Wf/H 


1.0 

3.432 

2.0 

3.996 

3.0 

4.254 

4.0 

4.401 

5.0 

4.497 

6.0 

4.565 

7.0 

4.615 

8.0 

4.654 

9.0 

4.685 

10.0 

4.710 

11.0 

4.770 

12.0 

4.789 

Table  (1Q.    Influence  of  1/H   on  predictions  based  on  equations  (11)  or 
.(16).     b  =  .0593,    7^=200,     D=  470  sec"1.     H  =  .  0643,  =417.16. 
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lt-» 

H-in 

Vo-ft/sec 

Wf/H 

0.  1 

.0643 

134.08 

323. 

055 

3.  325 

0.2 

ii 

3.698 

0.3 

ti 

3.852 

0.4 

it 

3.935 

0.5 

it 

3.988 

0.6 

ii 

4.025 

0.8 

ii 

4.071 

1.0 

ii 

4.  100 

0.  1 

.0644 

152.60 

417. 

16 

3.807 

0.2 

ii 

4.273 

0.3 

ii 

4.469 

0.5 

it 

4.646 

0.7 

it 

4.728 

1.0 

it 

4.793 

0.  1 

.0638 

180.06 

591. 

78 

4.534 

0.2 

ti 

5.  150 

0.3 

it 

5.416 

0.5 

it 

5.662 

0.  1 

.0638 

165.06 

497. 

30 

4.  173 

0.2 

H 

4.707 

0.3 

ti 

4.935 

0.5 

ii 

5.142 

0.7 

it 

5.239 

1.0 

ii 

5.316 

0.  1 

.0647 

233.  0 

963. 

55 

5.  594 

0.2 

ii 

6.499 

0.3 

ii 

6.912 

0.5 

ti 

7.303 

0.7 

it 

7.492 

^able  'IT). 

1.0 

it 
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7.644 

(continued) 

Table  (ID.    Influence  of  hinge  width  on  the  predictions  based  on 
equations  (11)  and  (16).     b  =  0.593,     D  -  500  sec"    ,    T^=  200. 
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APPENDIX  A 

i 

COMPUTER  PROGRAMME  FOR  TABLES 
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